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Abstract 

The coherence of two, coupled two-level systems, representing vibrational modes in a semiclassical 
model, is calculated in weak and strong fields for various coupling schemes and for different relative 
phases between initial state amplitudes. A relative phase equal to tt projects the system into a 
dark state. The selective excitation of one of the two, two-level systems is studied as a function of 
coupling strength and initial phases. 



I. INTRODUCTION 



Investigations of coherent dynamics in multidimensional systems induced bv ultrafast 
shaped pulses is being carried out in a number of research groups 

Special attention is paid to the development of femtosecond laser techniques for control 
over molecular motion in stimulated Raman spectroscopy, S, , Q, [ll | • In these works 
the selective Raman excitations are achieved with femtosecond laser pulses with spectral or 



phase modulation. The goal is to prepare a specific quantum state to facilitate 'unusual' 
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structural changes or chemical reactions, p, |9|, |12|, |13( . Another fundamental application is 
related to the development of quantum memory for information storage, A composition 
of two-level systems, e.g., vibrational normal modes within an ensemble of molecules, driven 
by an optical field, may serve as a storage device of quantum information. 

One of the important steps needed for efficient control of molecular motion is an under- 
standing of the factors that govern the system's time evolution. In this paper femtosecond 
pulse shaping is discussed that allow for the selective excitation of unresolved, coupled Ra- 
man transitions. 

In P|, a comparative analysis of Raman spectra of liquid methanol and a mixture of 
benzene and deuterated benzene showed experimental evidence for the dependence of the 
selective excitation on intramolecular coupling of normal vibrational modes in a molecule. 
In this paper a semiclassical model is developed for the interaction of a shaped ultrafast laser 
pulse with two, coupled two-level systems, representing coupled vibrational modes in a single 
molecule. Specific questions about the mechanisms of interaction of an external field with 
molecular media are addressed such as how the coupling via an external field influences the 
controllability of selective excitation and how the result depends on the coupling strength. 
We also investigate the ways of implementing a coupling mechanism for coherent control. 
Within our model selective, high level coherence can be achieved in two coupled two-level 
systems by a choice of the relative phase of the initially populated states. For a particular 
phase, a dark state is formed with zero eigenvalue of the interaction Hamiltonian. It is known 
that for a quantum system in a dark state the prevention of decoherence may be achieved, 
e.g., jl5| . Molecules prepared in such a state may be useful for quantum computation and 
communication. 
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II. BASIC FORMALISM 



A semiclassical model is used to describe the interaction of an ultrafast laser pulse with a 
molecular medium using stimulated Raman scattering. The model is represented schemati- 
cally in Fig.l, where two, two-level systems describe two Raman active modes in a molecule. 
Levels |1 > and |3 > are at zero relative energy, while levels |2 > and |4 > have energies 
T1UJ2 and hu4, respectively. Transition dipole moment matrix elements of the levels with a 
virtual intermediate state \b > are equal to fiib- Generally they may be different. We inves- 
tigate the effects in weak and strong fields caused by coupling between normal vibrational 
modes in a molecule and by the relative phase of the amplitudes of the initially populated 
states. Transition matrix elements of the 3-4 two-level system are assumed to be equal, 
f^3b = f^4b, and transition matrix elements of the 1-2 two- level system satisfy the following 
conditions — = — = r. The equations of motion for the probability amplitudes of two 
coupled two-level systems are: 
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where x is a time-dependent variable equal to I(t), I(t) is the pulse intensity envelope 
and A is the detuning of the frequency of the pulse from the frequency of the virtual state 
\b >. Note, that in our model the pulse intensity envelope I(t) is the same for all transitions. 
The specific form for the pulse shape is taken such that in the weak field regime the pulse 
selectively excites transitions of predetermined frequencies, while in strong fields the result 
depends on the intensity of the field 1$. The intensity envelope I(t) is defined as a real part 
of the Fourier transform of a function f{uj) specified as 



— e 



(1 



— e 



-(u;-W2) 2 7V 



(2) 



where a, T and T\ are free parameters. In the vicinity of resonances this function has the 
spectral profile identical to the one suggested in [16| except for the fact that it contains a 
"dc" component centered at u = 0, to ensure that I(t) is positive. 

The real part of the inverse Fourier transform of the spectral density (2) gives a temporal 
pulse function: 



I{t) = I C (l - A/Ccos{u 4 t) - B/Ccos{{uo 2 - Au^)t)j , (3) 

t 2 = T 2 + T?,Auj = uj 4 -uj 2 . 

The temporal profile of the pulse function is shown in Fig. 2 (a) for parameters corresponding 
to experimental data on the molecular gas C0 2 5], T = 3, Ti = 3, and u 2 — l,k>4 = 1.1, 
where are in frequency units, and T, 7\ are in inverse frequency units. The chosen 
parameters T, T\ give a pulse duration corresponding to an impulsive regime of interaction 
f|. In Fig. 2 (b) the Fourier transform I (to) of Eq.(3) is presented showing significantly 
larger intensity at frequency uj 4 than at uj 2 . In weak fields the solution of Eq.(l) is 

(4) 

a 2 = IZo I(t)e^dt = i^lM. 

As a consequence, when I(t) is used in Eq.(l), the solution in a perturbative limit leads to 
a weak excitation of coherence | P12 1 , and an efficient excitation of coherence \p3a\ with the 
magnitude proportional to I ((J 4). 

In strong fields the exact numerical solution of Eq.(l) shows that a choice of the field 
strength parameter Iq provides control over excitations resulting in maximum coherence for 
either |p 34 | or \p\ 2 \. These results are in agreement with those published in for two 
uncoupled two-level systems. 



III. MODE COUPLING 



A numerical solution of the time-dependent Schrodinger equation (1) describing two cou- 
pled two-level systems was obtained with I(t) given by Eq.(3) and two values of r, equal to 
I and 1. The results reveal the importance of the relative phase between the initial state 
amplitudes ai and a^. This relative phase could be established by optical pumping into 
state |1 > and using a Raman pulse to create the |1 > — 13 > state coherence. Coherences 
|pi 2 | and |p 34 | are calculated as an average over relative phases between initially populated 
states |1 > and |3 >. For r = 1, |p 12 | and |p 34 | are shown in Fig.3(a), (bold dashed and 
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bold solid lines, respectively), as a function of the dimensionless intensity of the ultrafast 
laser pulse. This phase averaging is equivalent to calculating |p 34 | and \p i2 \ as a sum of 
two contributions resulting from two initial population distributions: (pn = |, p 33 = 0) 
and (p 33 = \-, Pu = 0). This approach eliminates the role of the phase between initially 
populated states. Coherences as a function of the intensity of the field calculated in such a 
way are identical to those in Fig. 3 (a). Also shown in Fig. 3 (a), are values of |p 34 | and IP12I, 
(represented by thin solid and dashed lines, respectively), when there is no coupling between 
the two, two- level systems (obtained formally by setting r 2 = l,r = in Eq. (1)). From a 
comparison of the two sets of curves it is seen that phase averaged solution for two coupled 
two-level systems gives much lower values of | P12 1 and |p 3 4| than that for uncoupled systems. 
In weak fields the coherence |p 3 4| prevails over \p\2\ which is in agreement with Eq.(4). The 
coherence |pi 2 | increases faster than that for the case of two uncoupled two-level systems. 
In strong fields |pi 2 | and |p 34 | oscillate synchronously due to nonadiabatic interactions with 
the dc component of the field. The coherence |pi2| almost always is larger than |p 34 | owing 
to its slightly smaller transition frequency. 

Numerical results show that I(t) may result in significant differences between | 1 and 
|p 34 | in a system with different coupling constants p«. For r — \ the dependence of |p 34 | and 
|pi 2 | on the intensity of the field is shown in Fig. 3(b) by bold solid and bold dashed lines, 
respectively. For such coupling constants, the probability of population transfer between 
two-level systems is equal to 1/2, between states |1 > and |2 > it is equal to 1/4, and 
between states |3 > and |4 > to 1. As the result, population flows from the 1-2 to the 3-4 
two- level system, maximizing coherence |p 34 |. 

For various values of intensity of the field Iq the dependence of |pi2| and |p 34 | on the 
relative phase between initially populated states |1 > and |3 > has been considered. In 
Fig. 4(a) the case for I = 2.6257T, r = 1 is represented (which gives almost equal values of 
| yOi2 1 and |p 34 | for the phase averaged solution). For phases from zero to ir, |p 34 |/|pi2| > 1, and 
for phases from n to 2n, |p 34 |/|pi2| < 1- Consequently, using phase control of the initially 
populated states allows one to enhance the coherence between desired vibrational levels. 
The behavior of |p 34 |/|pi2| is sensitive to the intensity of field. For example, when J = 0.57T, 
(which gives a phase averaged solution |p 12 | = 0.163 and |p 34 | = 0.059 in Fig.3(a)), the phase 
dependent calculation shows that | P12 1 / 1 P34 1 > 1 f° r all phases except for <p = ir, see Fig. 4(b). 
The maximum value of coherence |p^ a:r | = 0.4 at = 7ix/A is much higher than that for 
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the phase averaged solution. When the initial relative phase is equal to 7r, coherences 
I/012 1 and I/034 1 are equal to zero. This is the case for any external field. Populations of all 
states exhibit no time evolution. This result indicates an existence of a dark state as can 
be deduced directly from the Hamiltonian in Eq.(l). For an arbitrary value r, the necessary 
conditions for a dark state are — = r 2 and 6 — ir. 

IV. SUMMARY 

We presented a semiclassical model of the coherent control of excitation of Raman transi- 
tions in two coupled two-level systems using a broad-band shaped laser pulse. We analyzed 
the effects caused by the coupling between four levels via a laser field. The coupling is 
shown to cause an efficient population transfer between two-level systems. Effects caused by 
the relative phase between the initially populated states were analyzed. Depending on the 
intensity of the field, the initial phase may cause predominance of the coherence of one two- 
level system with respect to another. When the relative phase between initially populated 
levels |1 > and |3 > is equal to n, two, two-level systems encounter stationary conditions 
regardless of the external field strength, implying the existence of a molecular dark state. 
These results may be useful for investigation of decoherence processes caused by the driving 
field. 
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FIG. 1: Schematic picture of a model system consisting of two, two-level systems having frequencies 
u>2 and CJ4. Initially, the lower levels are populated evenly. The transitions between four levels are 
driven by an off-resonant femtosecond pulse. 
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FIG. 2: (a) Intensity temporal profile I(t) for T = 3, T\ = 3 [uj 2 i}> (b) Fourier transform of I(t). 
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FIG. 3: Intensity dependence of the coherences of two, two-level systems; \p\2\ is shown by dashed 
lines and \pza\ by solid lines. In (a) bold curves depict the case for r = 1 corresponding to two 
coupled two- level systems with equal coupling constants; thin curves depict coherences for two 
independent two-level systems. The phase averaged solution gives much lower values of coherences 
than that for zero phase and coupling. In (b) bold curves show coherences for r = 1/2, and thin 
curves for r = 1. Weak coupling constants of the 1-2 system result in efficient population flow 
toward the 3-4 system, strongly coupled to the field. 
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FIG. 4: Coherence of the 1-2 and 3-4 two-level systems as a function of initial relative phase 
between levels |1 > and |3 > for r = 1 and equal initial populations of these levels, (a) Jo = 2.625-7T, 
(b) I = 0.5vr. 
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